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Background: Recent observational data on neutron star masses and radii provide stringent constraints on the 
equation of state of neutron rich matter [l] . 

Purpose: We aim to develop a nuclear energy density functional that can be simultaneously applied to finite 
nuclei and neutron stars. 

Methods: We use the self-consistent nuclear density functional theory (DFT) with Skyrme energy density func- 
tionals and covariance analysis to assess correlations between observables for finite nuclei and neutron stars. In 
a first step two energy functionals - a high density energy functional giving reasonable neutron properties, and a 
low density functional fitted to nuclear properties - are matched. In a second step, we optimize a new functional 
using exactly the same protocol as in earlier studies pertaining to nuclei but now including neutron star data. 
This allows direct comparisons of performance of the new functional relative to the standard one. 

Results: The new functional TOV-min yields results for nuclear bulk properties (energy, r.m.s. radius, diffraction 
radius, surface thickness) that are of the same quality as those obtained with the established Skyrme functionals, 
including SV-min. When comparing SV-min and TOV-min, isoscalar nuclear matter indicators vary slightly while 
isovector properties are changed considerably. We discuss neutron skins, dipole polarizability, separation energies 
of the heaviest elements, and proton and neutron drip lines. We confirm a correlation between the neutron skin 
of 208 Pb and the neutron star radius. 

Conclusions: We demonstrate that standard energy density functionals optimized to nuclear data do not carry 
information on the expected maximum neutron star mass, and that predictions can only be made within an 
extremely broad uncertainty band. For atomic nuclei, the new functional TOV-min performs at least as well as 
the standard nuclear functionals, but it also reproduces expected neutron star data within assumed error bands. 
This functional is expected to yield more reliable predictions in the region of very neutron-rich heavy nuclei. 



PACS numbers: 21.60.Jz,21.65.Cd,26.60.-c,26.60.Kp 

I. INTRODUCTION 

We encounter the same material both in the labora- 
tory and in outer space. This basic principle has been 
remarkably powerful for centuries. For Newton, the ma- 
terial was mass, and he observed that gravity acted in 
the same way on an apple and on the moon. In the 19th 
century, the materials were the chemical elements. As- 
trophysics grew from spectral lines being the same in the 
laboratory and in the cosmic space. Indeed, the second 
element helium is named after the Greek God of the sun. 
The material, we focus on in this paper, is neutron rich 
matter. Both in the laboratory and in astrophysics, it 
has the same neutrons, the same strong interactions, and 
the same equation of state (EOS). We aim to develop a 
universal nuclear energy functional that can be applied to 
stable, neutron rich, and very heavy nuclei, and to neu- 
tron rich matter in astrophysics including neutron stars. 

A frequently adopted strategy is to optimize an energy 
density functional (EDF) to properties of finite nuclei. 
However, this can leave many properties of bulk neutron 
matter poorly constrained. Alternative optimization pro- 
tocols are based on considering both experimental data 



on finite nuclei and theoretical pseudo-data on nuclear 
and neutron matter. For example, the SLy4 [5] EDF was 
fit to both nuclei and the variational Fermi hypernetted 
chain calculations of neutron matter [SJ. Chamel et al. 
calculated the maximum mass of neutron stars for three 
EOS fit to finite nuclei and microscopic neutron matter 
calculations [4], and such analysis helped them to so se- 
lect one EDF (BSk21) that does well on both kinds of 
data. However, it is to be noted that current microscopic 
calculations of neutron matter are highly model depen- 
dent; in particular, the impact of the poorly known three- 
neutron forces (especially their T = 3/2 component) 
seems to be crucial [5J [B]. Furthermore, four-nucleon 
forces [7j may also be important, especially at high den- 
sity. Indeed, at this time, our ability to calculate neu- 
tron matter properties at high densities is fundamentally 
limited. For example the chiral effective field theory cal- 
culations of Ref. [8] do not converge at densities much 
beyond normal nuclear (saturation) density. Although 
microscopic approaches with phenomenological two- and 
three-nucleon interactions such as those of Ref. [5] do not 
explicitly depend on a chiral expansion, it is not clear 
how can they avoid large ambiguities that are present in 
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a chiral approach at densities beyond where the chiral 
expansion converges. 

Recently, there has been considerable progress in de- 
termining neutron star (NS) masses and radii [lj. Dc- 
morest et al. have accurately measured a 1.97 ±0.04 M Q 
neutron star using Shapiro delay of the radio signal [3]. 
This model-independent result immediately provides an 
important lower limit for the maximum mass of a neu- 
tron star. The EOS of neutron rich matter must have a 
high enough pressure to support at least this mass against 
collapse into a black hole. In addition, neutron star radii 
and masses have been inferred from X-ray observations 
of both quiescent stars in globular clusters and stars un- 
dergoing thermonuclear bursts [TDJ [TT] . These results de- 
pend on atmosphere models and on the model assumed 
for the X-ray bursts |12j . For example, Suleimanov et 
al. employ more sophisticated atmosphere models [T3] 
and obtain a larger radius than Steiner et al. lTj. while 
Ref. 1 14] finds that extracted neutron star radii depend 
on the assumed composition of the atmosphere. In the 
future, one may be able to obtain neutron star radii from 
observations of X-ray burst oscillations [TS] ■ 

In principle heavy ion collisions can provide data on 
the symmetry energy at high densities and help constrain 
energy functionals [IB]. For example, Ref. [T7] studied 
EOS by means of kaon production at subthreshold en- 
ergies. However, the interpretation of present heavy ion 
measurements may be model dependent. Finally, antici- 
pated gravitational wave observations |18j from neutron 
star mergers may provide constraints on EOS [T5] . 

Modern energy functionals such as UNEDFO [20 , UN- 
EDF1 [21 , and SV-min [22] have been optimized to a 
large variety of nuclear data. (For reviews of Skyrme 
Hartree-Fock (SHF) calculations of nuclei see for exam- 
ple Refs. [231121].) Ref. [25] provides an example of an 
application of SHF approach to neutron star structure 
while Ref. [26] also reviews Skyrme EDFs and discuss 
their divergent predictions for neutron and nuclear mat- 
ter. (Note that some of the constraints imposed in their 
review may be model dependent.) 

Agrawal et al. |27j optimized a relativistic mean field 
(RMF) Lagrangian to nuclei and some assumed neutron 
matter properties. Fattoyev et al. [28\ discussed opti- 
mization of isovector parameters of both RMF and SHF 
models using microscopic calculations of neutron matter. 
Other RMF functionals have also been determined while 
paying attention to some NS properties [29] [30]. Low 
density matter at sub-saturation density has been consid- 
ered in Ref. [3T] . Finally, RMF calculations at high den- 
sity, coupled with virial expansion and nuclear statistical 
equilibrium calculations at low density, have been used 
to generate astrophysical EOS [221 [33]; these give the 
pressure as function of density, temperature, and proton 
fraction, and can be used in simulations of supernovae, 
neutron star mergers, and black hole formation. 

We implicitly assume that neutron rich matter can be 
at least approximately described by a single EOS from 
nuclear density to the central density of massive neu- 



tron stars. We do not exclude phase transitions to non- 
nucleonic degrees of freedom (exotic components), but 
consider an EOS that changes smoothly with density. 
For instance, a high density phase with hyperons has 
been considered in Refs. [331 US], and a hadron-quark 
crossover was considered in Refs. [56] 137]. In general the 
appearance of exotic components could reduce the pres- 
sure of neutron rich matter. However, the observation of 
a massive 1.97 M Q neutron star [5J strongly suggests that 
any softening of the EOS at high densities associated with 
phase transitions is modest. Otherwise the EOS might 
not be able to support 1.97 M against collapse to a 
black hole. Note that quark matter may have an EOS 
similar to that of pure hadronic matter [38]. In this case, 
NS mass and radius observations could not distinguish 
a quark-hadron hybrid star from a purely hadronic star. 
In either case, the assumption of a single EDF should be 
applicable. 

In this work, we choose to optimize the nuclear EDF 
to both nuclear data and neutron star masses and radii. 
This has the advantage of treating finite nuclei and neu- 
tron star observables on an equal footing. We consider 
a mixed set of fit-observables consisting of varios bulk 
properties of spherical nuclei and NS data, including the 
maximum mass and the radius of a 1.4 Mq star. A new 
functional is obtained by optimizing the coupling con- 
stants to this extended dataset. This approach allows 
one to see how the NS data impact various parts of the 
energy functional and predictions in the neutron-rich ter- 
ritory, especially around the neutron drip line. 

Our paper is organized as follows. Section [Tl] overviews 
our model. In Sec. |II A| we describe our original SV-min 
functional that has been optimized to nuclear properties 
only. Next in Sec. |II B| we review the TOV equations 
of general relativity that describe NS structure. In Sec. 
|II C| we describe an intermediate procedure that simply 
matches a high density energy functional, that gives rea- 
sonable NS properties, to a low density functional ad- 
justed to nuclear properties. In Sec. IIP we describe 
how we optimize EDF to both finite nuclei and NS prop- 
erties simultaneously. Results of this optimization are 
presented in Sec. Ill The NS observables depend on 



properties of very neutron rich matter. This should al- 
low our functional fit to NS to make more accurate pre- 
dictions for very neutron rich and for very heavy nuclei. 



Therefore, in Sec. IV we apply our functional to nuclei 



near the drip lines over a very large range of mass num- 
ber A up to superheavy systems. Finally, Sec. [Vjcontains 
the conclusions of our work. 



II. ENERGY DENSITY FUNCTIONAL 
FORMALISM 

A. Existing functionals fit to finite nuclei 

This section briefly outlines the nuclear DFT in the 
self-consistent SHF variant, and the optimization strat- 
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egy employed in this study. (For an in-depth presenta- 
tion, we refer the reader to Refs. [23] [23] and references 
cited therein) . The main ingredient of the SHF theory is 
an expression for the total energy: 



E = E, 



kin 



J d 3 f £sk + ^Coul + ^pair — 



(i) 



where -Ekin is the kinetic energy; £sk ~~ the Skyrme EDF; 
EqouI ~ the Coulomb term (where the direct Hartree en- 
ergy is treated exactly and the exchange term in Slater 
approximation) ; E pa i r is the pairing functional; and -E CO rr 
is the correlation term that accounts for dynamical corre- 
lations (cm. correction, rotational correction). The key 
piece in the SHF approach is the Skyrme EDF, which 
represents the effective nuclear interaction between nu- 
cleons and can be partly derived by considering a low 
momentum expansion of the density matrix [39MT] : 



£sk= ]T C p T {l + D p T p a ) p 2 T + c£ p PT A PT 



T=0,1 



PT T T + Crp^rp + Crp PT^'^Ti (2) 



where 



Pt=o = Pp + Pn 



Pt= 



Pp - Pn 



(3) 



are the isoscalar and isovector densities, respectively. 
The local densities in the functional are the particle den- 
sity p q , kinetic-energy density T q and spin-orbit density 
J g , where the isospin index q = p,n labels protons and 
neutrons (for details and definitions of the coupling con- 
stants, see [201 H31 EH)- The functional Q shows, in 
fact, only the time-even couplings which suffice for cal- 
culations of ground states (g.s.) of even-even nuclei. The 
form ([2]) is used for the parameterizations SV-min [32] 
and TOV-min, which will be introduced here. We will 
also consider an alternative to the standard density de- 
pendence (oc p 2+a ) in terms of Pade approximants yield- 
ing the parameterization RD-min [43] . 

Lacking still sufficiently precise input from ab-initio 
many-body theories, we adjust the parameters of the 
Skyrme functional phenomenologically. To that end, we 
have scrutinized nuclei for correlation effects [33] and cho- 
sen a large sample of spherical nuclei which have negligi- 
ble correlations thus being well described in pure mean- 
field theory. The compilation of fit-observables used for 
the optimization of SV-min and TOV-min is given in Ta- 
ble m 

The model parameters are adjusted to these data by 
virtue of a least-squares fit, for details see Ref. [32] • The 
unconstrained fit of the functional |2]) to the dataset from 
Table [I] yields the parameterization SV-min. Note that 
only information from finite nuclei enters its calibration. 
A similar fit using the same data but a modified density 
dependence yielded EDF RD-min [33]. 

The well developed technique of least-squares covari- 
ance analysis yields not only optimized parameterizations 
which can be used to extrapolate to unknown regions, but 
it also allows to estimate an extrapolation error. This 



TABLE I. Compilation of phenomenological input for the op- 
timization of SV-min 22 and TOV-min EDFs. Abbreviations 
used are: BE - g.s. binding energy; i?diff _ charge diffraction 
radius; a - charge surface thickness; r c h - charge r.m.s. ra- 
dius; e; s - spin-orbit splitting of selected single-particle states; 
NS - NS data (only for TOV-min). 
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is illustrated in Fig. [T] which shows the binding energy 
curves in homogeneous nuclear matter and neutron mat- 
ter for SV-min, RD-min, and TOV-min together with 
their extrapolation uncertainties. Note that the uncer- 
tainty is increasing dramatically with increasing density 
for SV-min, which was calibrated to finite nuclei only (cf. 
Ref. [45] for more discussion). As the neutron matter 
EOS is the main ingredient for the description of neu- 
tron stars, predictions for neutron stars with standard 
Skyrme EDFs (here SV-min and RD-min) are bound to 
be plagued by large uncertainties, see discussion below. 



B. TOV equations and neutron star structure 

The structure of a non-rotating cold neutron star can 
be characterized by the relationship between the gravi- 
tational mass M and the radius R of the star [TJ. This 
relation is obtained by solving the Tolman-Oppenheimer- 
Volkov (TOV) equations [461E7]: 



dp(r) 
dr 



Ge(r)M(r) 
c 2 r 2 

4irr 3 p(r) 
1+ c 2 M(r) 



1 



P(r) 
e(r) 

2GM(r) 



(4) 



M(r)=4n r /2 dr' >(/) = 4tt / r' 2 dr' e(r' ) / c 2 . (5) 
Jo Jo 

Here p{r) — e/c 2 is the mass density at the distance r and 
e(r) is the corresponding energy density. The enclosed 
mass M(r) is the mass inside a sphere of radius r, and p 
is the pressure. 
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FIG. 1. (Color online) Binding energy per nucleon in ho- 
mogeneous neutron matter (a) and symmetric matter (b) as 
a function of density. Shown are results for the parameteri- 
zations SV-min [22], RD-min [43], and TOV-min (this work). 
For SV-min and TOV-min, we show also error bars represent- 
ing the uncertainty in the extrapolation. The error bars for 
RD-min are comparable to those of SV-min and dropped for 
better readability of the figure. 



The main physics ingredient necessary to solve equa- 
tions Q and ([5]) is the equation of state (EOS) providing 
the pressure as a function of the energy density p = p(e) . 
The coupled equations for p(r) and M{r) can be inte- 
grated starting from r — 0, and a value for the pressure at 
the center of the star po, up to the point where p{R) = 0. 
The result is the mass radius relation M — M(R), which 
is expected to be consistent with our assumptions of a 
maximal mass of 2.2Af and a radius of 12.5km for a 
1.4M neutron star, see below. As a neutron star is 
not only made of neutrons, but is a mixture of neutrons, 
protons, electrons, and possibly muons - due the weak 
decay of neutrons and electron capture processes on pro- 
tons (for higher pressure even hyperons or kaons) - the 
EOS has to be calculated for /3-equilibrium [35] • Neutron 
star matter is characterized by the following processes: 



P 



p + fi 



(6) 



The corresponding chemical potentials should fulfill the 
/3-equilibrium conditions: 

Mn = Mp+Me) Mm = Me- (7) 

Here each chemical potential is given by the canonical 



de 



(8) 



with the total energy density e containing all nucleon and 
lepton contributions for the particle number densities rij 
(for details see [H]). In addition, charge neutrality has 
to be fulfilled: 

n p =n e + V (9) 

Lepton chemical potentials are calculated for nonintcr- 
acting Fermi gases. In the following, all EOS are calcu- 
lated at /3-equilibrium. 
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FIG. 2. (Color online) Mass-radius relation of SLy4 [2], 
UNEDF0 [20] and SV-min [25]. The uncertainty band for 
SV-min is shown. This band is estimated by calculating the 
covariance ellipsoid for the mass M and the radius R at each 
point of the SV-min curve as indicated by the ellipsoid. Also 
depicted (dotted lines) are uncertainty limits for TOV-min. 

The results for the mass-radius relation of neutron 
stars M(R) is shown in Fig. [2]for SV-min, SLy4, and UN- 
EDF0. In addition, the estimated uncertainty band for a 
prediction using SV-min is shown. As both observables 
are obviously correlated it is not possible to estimate the 
uncertainty for mass M and radius R separately. For this 
reason, the error band is obtained by calculating the co- 
variance ellipsoid for mass M and radius R of a neutron 
star for each point of the M(R) curve as indicated by 
the ellipsoid in Fig. [2] This ellipsoid is a measure for the 
correlation of two observables as presented in [23], and 
provides information about the uncertainty of an observ- 
able in dependence of a second one. The area covered by 
all ellipsoids can be interpreted as the error band for a 
prediction using SV-min. 

We now discuss reasonable expected values for the 
maximum neutron star mass M max and the radius of a 
1.4M neutron star R\a- Demorest et dl. have observed 
a 1.97 ± 0.04 M Q star [9], so clearly M max > 1.97M . 
However elementary population synthesis considerations 
imply that the observation of a 1.97M star suggest that 
the actual maximum mass must be, at least, somewhat 
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greater than this value. This will allow for a reasonable 
probability to observe a 1.97-Mq star, if one draws neu- 
tron star masses from a realistic distribution. Therefore 
we assume that 



Mr, 



2.2 ± 0.2 M, 



■ 



(10) 



This choice for M max also provides room to accommodate 
the observation of somewhat more massive stars in the 
future. 

X-ray observations can provide information on neu- 
tron star radii. For example, from luminosity and sur- 
face temperature measurements one can infer an emit- 
ting area. However, the interpretation of X-ray obser- 
vations may be model dependent. One can be sensitive 
to the employed neutron star atmosphere model and or 
the assumed model for photosphere radius expansion X- 
ray bursts. The radius of a 1.4 Mq star R\a between 
10.4 and 12.9 km has been inferred in Ref. [H]. How- 
ever Refs. [13] use more sophisticated neutron star at- 
mospheres to model a long X-ray burst and infer larger 
radii. At this point, we believe the final word on X- 
ray observations of neutron star radii has not yet been 
written. Future observations and/or more sophisticated 
theoretical interpretations may change the inferred radii. 
Therefore, at this time, we somewhat arbitrarily adopt a 
value for R±a near the upper end of range of Ref. [48] : 



R 1A = 12.5 ±0.5 km. 



(11) 



Here, the assumed 0.5 km error is somewhat arbitrary. 
This value insures that fits to both nuclear and NS data 
give reasonable weight to the neutron star observations, 
see below. 

As shown in Fig. [2] the established Skyrme functionals 
can not reproduce the expected M max ( 10 ) and R\a (11 1- 



However, these values are at the edges of the error band, 
and they can be reached by a new EDF fit that includes 
NS data. Note that we do not explicitly include an EOS 
for the nonuniform neutron star crust. Instead we sim- 
ply use the EOS for uniform matter implied by a given 
Skyrme interaction even at low density. This error has 
essentially no impact on M max . However, it somewhat 
reduces the radius of 1.4M stars compared to predic- 
tions that employ realistic crust EOS at low densities. 
Furthermore, this approximation significantly underesti- 
mates the radius of very low mass stars. 



C. Matching to a high density equation of state 

As a first step in creating a Skyrme functional for pre- 
dicting properties of both finite nuclei and NS data, we 
match an established nuclear Skyrme functional to an 
ansatz for the equation of state at high densities, for 
which we take LS220 [49]. To deal with neutron stars, it 
is necessary to consider a wide range of densities, starting 
from below the nuclear saturation density pq « 0.16 fm -3 
up to as high as p ~ 10 po- The question is what is the 



optimal density to switch between the two equations of 
state. To answer this question, the following ansatz for 
a combined equation of state can be made 

Pc(p) = /(p)Pedf(p) + (1 - f(p))P hi M, (12) 



"Pc(p'), , 

dp , 



where a switch factor 



p-ipo 

f(p)=[l + e •/> 



(13) 



(14) 



ensures a smooth matching. Here -Pedf(p) is the pres- 
sure for a Skyrme energy functional and -Phigh(p) is the 
pressure for the assumed high density functional. In the 
switch factor, a p — 0.2 p is a fixed diffuseness parame- 
ter and £ is a parameter that determines the matching 
density. The advantage of this matching procedure is 
that any energy functional can be used to describe fi- 
nite nuclei, while still providing a reasonable description 
of neutron stars. Figure [3] shows the neutron star mass 
versus radius relation from combining the UNEDF0 and 
LS220 functionals, with several values of £. The maxi- 
mum mass is seen to be relatively insensitive to £. Based 
on this result one can conclude that the low-density, or 
nuclear, part of the EOS carries no information on Af max . 
That is, it makes little sense to scrutinize existing nuclear 
EDFs with respect to this quantity. On the other hand, 
the radius of a 1.4M neutron star can vary by up to 
1 km depending on the transition density used to match 
the two equations of state. 
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FIG. 3. (Color online) The mass-radius relation of an EOS 
matching the EDF UNEDF0 at low densities with the func- 
tional LS220 at high densities for matching densities p = £po 
with £ = 2.0, 2.5, 3.0, and 3.5 (fill 
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D. Optimizing EDF to nuclear and NS data 

In this section we describe optimization of a single EDF 
to both nuclear data listed in Table U and to the maxi- 
mum mass M max and radius R\a of neutron stars, as in- 
dicated in Fig. [2j We use exactly the same optimization 
protocol as in Rcf. [22 . As a starting point, we choose 
the Skyrme EDF SV-min. To include the additional in- 
formation for neutron stars, the TOV equations Q and 
([5]) are solved at each step of the optimization procedure, 
assuming matter in /3-equilibrium. The result is the new 
Skyrme EDF TOV-min (see Table [HI). Figure [4] presents 



TABLE II. Optimized parameter set of TOV-min expressed 
in terms of the traditional (t, x) parameterization of the 
Skyrme force and in terms of C coupling constants. 



to (MeV-fm 3 ) 


-2129.735 


CZ (MeV-fm 3 ) 


-798.650 


ti (MeV-fm 5 ) 


305.398 


(MeV-fm 3 ) 


356.703 


ta (MeV-fm 5 ) 


362.532 


C PP D (MeVfm 3+3a ) 


872.316 


t 3 (MeV-fm 3+3a ) 


13957.064 


C PP D (MeVfm 3+3t *) 


-524.922 


x 


0.169949 


C pT (MeVfm 5 ) 


9.02926 


Xl 


-3.399480 


C pT (MeV-fm 5 ) 


52.5831 


X2 


-1.782177 


CS Ap (MeV-fm 5 ) 


-55.0048 


%3 


0.402634 


C pAp (MeV-fm 5 ) 


-97.5411 


64 (MeV-fm 5 ) 


37.07357 


C^ J (MeVfm 5 ) 


-79.1223 


b' 4 (MeV-fm 5 ) 


84.09737 


C pVJ (MeVfm 5 ) 


-42.0487 


a 




0.250388065 




Vp air ,p (MeVfm 6 ) 




630.516 




Vp air , n (MeVfm 6 ) 




629.072 




PO,pair (fm" 3 ) 




0.201322348 




& p (MeV-fm 2 ) 




20.7498207 




2^ (MeV-fm 2 ) 




20.7212601 





the mass-radius relation for TOV-min together with re- 
sults for SV-min and RD-min. The new functional per- 
forms well in reproducing the expected maximum mass 
and radius within the adopted errors (indicated by ar- 
rows). Note also that RD-min, using the modified den- 
sity dependence, also fits NS data without considering 
NS information during the optimization process. 



III. TOV-MIN PERFORMANCE 

In this section results for the TOV-min EDF are pre- 
sented and compared to results for some existing Skyrme 
functionals. 
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FIG. 4. (Color online) The mass radius relation for neutron 
stars obtained with SV-min, RD-min, and TOV-min. 



A. TOV-min properties 

To study the performance of the new functional, Table 
|III| lists the root mean square errors of SV-min, RD-min, 
and TOV-min on nuclear fit-observables: g.s. binding 
energy BE, charge diffraction radius -ffdiff , surface thick- 
ness a and r.m.s. charge radius r c h, and pairing gaps 
for protons and neutrons (A p and A„). Adding NS data 



TABLE III. Root mean square errors on nuclear fit- 
observables obtained by the parameterizations SV-min, RD- 
min, and the new functional TOV-min. 



r.m.s. error 


SV-min 


RD-min 


TOV-min 


BE (MeV) 


0.62 


0.63 


0.66 


-Rdiff (fm) 


0.028 


0.030 


0.028 


a (fm) 


0.022 


0.022 


0.022 


r ch (fm) 


0.014 


0.013 


0.014 


A p (MeV) 


0.11 


0.11 


0.11 


A„ (MeV) 


0.14 


0.14 


0.14 


Total x 2 


53.22 


55.90 


54.95 



has only a minor impact on the binding energy, while 
there is no change in performance for other nuclear bulk 
properties. This results in only a small increase in the 
total x 2 and shows that including NS data has almost no 
influence on nuclear bulk properties considered. 

Table |IV| shows nuclear matter properties as defined 
in Ref. [2UJ HI]- I n addition, the surface energy coeffi- 
cient a sur f and surface-symmetry energy coefficient a ssym 
are shown. These are calculated from a leptodermous ex- 
pansion of the energy functional (for details see Ref. 50 ) . 
Results for isoscalar properties (/9 eq , E/A, K, m* Jm and 
flsm-f ) vary slightly between SV-min and TOV-min, while 
results for isovector properties (a sym , K, L sym and a ssym ) 
differ considerably as a result of neutron star constraints. 
Of particular interest in the context of neutron struc- 
ture are: the symmetry energy £2 at saturation density, 
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TABLE IV. Nuclear matter properties as defined in [22] for 
SV-min, TOV-min, RD-min, and UNEDFO. For the isovector 
properties k, a sym , and L sym the calculated uncertainties are 
indicated. The surface energy coefficient a sur f and surface- 
symmetry energy coefficient a ssym are also given. 





SV-min 


TOV-min 


RD-min 


UNEDFO 


po (fm" 3 ) 


0.1610 


0.1610 


0.1611 


0.1605 


E/A (MeV) 


-15.91 


-15.93 


-16.11 


-16.06 


K (MeV) 


222 


222 


231 


230 


m* /m 


0.95 


0.94 


0.90 


0.90 


a sym (MeV) 


30.7±1.9 


32.3±1.3 


32.1±2.1 


30.5 


K 


0.08±0.29 


0.21±0.26 


0.04±0.32 


0.25 


L sym (MeV) 


45±26 


76±15 


60±32 


45 


a surf (MeV) 


17.6 


17.6 


17.6 


18.7 


a ssym (MeV) 


-51 


-44 


-55 


-44 



S2(po), and the slope of the symmetry energy 



3po 



dp 



(15) 



PO 



The values of a suxt and L sym in Table IV 



are consistent 



(within error bars) with the experimentally allowed range 

DEBET]. 



B. Correlations, neutron skin and dipole 
polarizability 



Figure [5] presents correlation between M max or Ria 
and nuclear matter properties predicted by SV-min. 
(Note that we show correlations for SV-min, rather than 
for TOV-min, because Ria and M max were explicitly 
constrained in optimizing TOV-min.) It is seen that Ria 
is very strongly correlated with a sym and L sym while the 
correlation is slightly weaker for M max . Other nuclear 
matter properties, in particular the saturation density 
po, correlate weakly with the NS properties. 

Table [V] displays the neutron skin thickness i? s kin = 
r™ s — r T p ms and electric dipole polarizability ao of 208 Pb 
as predicted by SV-min, RD-min, and TOV-min. Both 
observables are the focus of recent measurements and 
constrain the isovector part of the nuclear EDF that 
is known to be poorly constrained by previous data 
[151 [SH [55J [57]. For example, the PREX experiment at 
the Jefferson Laboratory has measured the neutron ra- 
dius of 208 Pb using parity-violating electron scattering 
[55HSJJ> while the electric dipole polarizability of 208 Pb 
was recently accurately determined to be — 14.0 ± 
0.4fm 2 /MeV at RCNP in Osaka in a high-resolutionn 
(Pip) measurement [52]. The predicted values of i? s ki n 
and cud in Table [V] are consistent with the data and the 
most recent theoretical estimates [53]. Since TOV-min 
yields larger value of L sym as compared to SV-min and 
RD-min, it is not surprising to see that it predicts an 
increased neutron skin thickness and cvd- Indeed, i? s kin 
and qd are strong isovector indicators that are well cor- 
related with i sym [451154] . 



A quantity characterizing the correlation between two 
observables A and B within a given model is the Pearson 
product- moment correlation coefficient Cab [52] . where 
a value Cab = 1 means the observables are fully corre- 
lated and Cab = implies they are totally uncorrelated. 
The covariance analysis in terms of Cab has been proven 
useful in previous surveys [551 SSI I53"H56] . 



TABLE V. Neutron skin i? s km and electric dipole polariz- 
ability «D of 208 Pb calculated using SV-min, RD-min and 
TOV-min. 





R skin (fm) 


q d (fm 2 /MeV) 


SV-min 


0.170±0.036 


14.2 ±0.7 


RD-min 


0.189±0.043 


14.1 ± 0.9 


TOV-min 


0.205±0.021 


15.0 ±0.5 
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FIG. 5. (Color online) Correlation Cab between selected 
nuclear matter properties and (a) the maximum mass of a 
neutron star and (b) the radius of a 1.4Mq neutron star ob- 
tained for the Skyrme functional SV-min. 



The covariance analysis allows one to study correla- 
tions between nuclear observables and NS properties. 
Figure [6] shows the covariance ellipsoid for i? s kin of 208 Pb 
and R\a calculated using SV-min. The corresponding 
correlation coefficient is large, Cab — 0.82, confirming 
a correlation between those two observables. An even 
stronger correlation, Cab = 0.95, was obtained in Ref. 
[55] using the relativistic EDF FSUGold. In general, a 
higher pressure for neutron matter near po increases both 
the neutron skin thickness, as neutrons are pushed out 
against surface tension, and Ria- However, Ria also de- 
pends on the pressure of neutron matter at high densities 
[63] . Therefore RD-min, which has a different density de- 
pendence, predicts almost the same Ria with a smaller 
-Rskin compared to TOV-min. 

Note that other neutron star properties have been pre- 
viously correlated with i? s ki n , see for example j56j l64] - 
The transition density from solid crust to liquid core was 
found to be anticorrelated with i? s ki n 56, 65 . In addi- 




0.8 0.9 1.0 1.1 1.2 



R(1.4 MgJ/R 

FIG. 6. (Color online) The covariance ellipsoid for -R s km 
in 208 Pb and the radius of a 1.4M NS calculated using SV- 
min. The mean values for SV-min are: Ria = 10.18 km and 
fiskin = 0.170 fm. Also shown is the correlation line corre- 
sponding to FSUGold [56] obtained using Ria = 12.66 ± 
0.46 km and R Bkin = 0.207 ± 0.037 fm. 

tion, the threshold density for the rapid cooling of neu- 
tron stars via the direct URCA process was found to 
decrease with i? s ki n |55] - 

IV. APPLICATIONS TO HEAVY NUCLEI AND 
DRIP LINES 



\ Fm 


SV-MIN — — - 
OV-MIN 
NEDF0 

exp o 







140 160 180 200 220 240 
neutron number 



FIG. 7. (Color online) Calculated (SV-min, UNEDF0, and 
TOV-min) and experimental two-neutron separation energies 
of even-even fermium isotopes. The error bars on SV-min 
and TOV-min results indicate statistical errors due to the 
uncertainty in determining their coupling constants. 

As seen in Sec. |IIIB| adding NS data helps constrain 
isovector nuclear matter properties and the properties 
of neutron-rich nuclei such as the neutron skin thick- 
ness. These additional data should allow one to make 
more accurate predictions for very neutron-rich nuclei. 



In this section, we study the effect of NS data on the 
structure of nuclei near the drip lines. Figure [7] shows 
two-neutron separation energies of even-even fermium 
isotopes using SV-min, UNEDF0, and TOV-min. For 
isotopes where experimental data is available, the mod- 
els agree and nicely reproduce the experimental data as 
discussed in |67j . Moving to more neutron rich isotopes 
the discrepancy between various predictions grows and 
ends up in different predictions for the location of the 
neutron drip line (S2 n =0). In the same way, statistical 
error bars for SV-min and TOV-min grow with increas- 
ing neutron number. However, TOV-min shows smaller 
error bars compared to SV-min because of the additional 
constraints on NS properties. Indeed the theoretical un- 
certainty band for the mass-radius relation is dramati- 
cally reduced when going from SV-min to TOV-min, see 
Fig.0 

The functionals SV-min and UNEDF0 show the phe- 
nomenon of re-entrant binding caused by shell effects |67j 
where some heavier isotopes are bound beyond the first 
two-neutron drip line. However the separation energies 
for TOV-min are slightly larger and do not exhibit the 
re-entrant behavior. 

A global survey of the two-nucleon drip lines predicted 
with TOV-min up to proton number Z = 120 is sum- 
marized in Fig. |U The TOV -min results are com- 
pared to the uncertainty band of a systematic study 
[57] (2012 Benchmark) using a variety of different nu- 
clear EDFs. Compared to the uncertainty band, the two- 
neutron drip line for TOV-min is slightly shifted towards 
more neutron-rich nuclei, and it lies near the outer bor- 
der of the error band of 2012 Benchmark. A reason for 
deviations of TOV-min from the error band around shell 
closures at N = 184 and N = 258 could be the non- 
existence of re-entrant binding for TOV-min as seen in 
Fig. [7] due to an increased binding (and 6*2,1) of neutron- 
rich nuclei. For the two-proton drip line, the TOV-min 
results agree very well with 2012 Benchmark. 



V. SUMMARY AND CONCLUSIONS 

In this study, we proposed and investigated new 
approaches for predictions of NS properties based on 
Skyrme functionals. We have shown that standard 
Skyrme functionals, adjusted to properties of finite nu- 
clei, can not reproduce the expected maximum neutron 
star mass M max and the radius of a 1.4M Q neutron star 
Ria, and that predictions can only be made within an 
extremely broad uncertainty band. 

In a first step, we calculated a new EOS by matching 
a high-density EDF, that gives reasonable NS properties, 
to a low-density EDF adjusted to nuclear properties. Us- 
ing this approach, by matching LS220 to UNEDF0, it 
is possible to predict a maximum mass of about 2.2M , 
but the radius of a 1.4M neutron star is relatively small, 
ranging from 10 to 11 km. Based on this exercise, we can 
conclude that the low-density part of the EOS carries no 
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FIG. 8. (Color online) Two-nucleon drip line of TOV-min (dark blue) compared to the uncertainty band from Ref. |67| 
(2012 Benchmark) obtained by averaging results of six different EDF: SkM*, SkP, SLy4, SV-min, UNEDFO, and UNEDF1. 
The dashed grey grid lines show the magic numbers known around the valley of stability (20, 28, 50, 82, 126) as well as the 
predicted regions of enhanced shell-stability in superheavy nuclei around N = 184 and 258. All 767 experimentally known 
even-even isotopes are shown with the stable nuclei indicated as black squares and the radioactive nuclei as green squares. 



information on M max . Therefore, EDFs optimized to nu- 
clear data cannot be used to predict this quantity, and 
scrutinizing existing functionals with respect to M mm 
makes little sense. 

In a next step, we have optimized a new Skyrme func- 
tional TOV-min by simultaneously considering nuclear 
and NS data. Results for standard nuclear bulk proper- 
ties are as good as results obtained with the established 
functionals. Nuclear matter properties vary slightly for 
isoscalar properties between SV-min and TOV-min while 
isovector properties are considerably changed. This is 
not surprising because the NS data mainly constrain the 
isovector channel. In particular, the values of a sur f and 
L sy m of TOV-min are increased with respect to SV-min, 
but they are still consistent with the experimentally al- 
lowed range [TJ [16] . Interestingly, the RD-min functional 
that has a modified density dependence does also very 
well on M max and -R1.4, in spite of the fact that it has 
not been optimized to NS data; its values of a sur f and 
L sym are fairly close to those of TOV-min. 

A quantity that contains information on the isovector 
channel is the neutron skin thickness. Adding NS data 
leads to an increased neutron skin for TOV-min com- 
pared to SV-min. Furthermore, a previously reported 
correlation between neutron skin and the NS radius is 
confirmed. We further checked the predictions of TOV- 
min for neutron rich nuclei and for the position of two- 
particle drip lines. Results for separation energies based 
on the new functional agree with established Skyrme 
functionals, such as SV-min or UNEDFO, in the regions 



where experimental data are available. The two-neutron 
drip line predicted by TOV-min is shifted towards the 
more neutron-rich side of the 2012 Benchmark uncer- 
tainty band [57] , while the two-proton drip line of TOV- 
min lies within the 2012 Benchmark. 

The new functional TOV-min is a first example of EDF 
optimization using both nuclear and neutron star data. 
While the results obtained with TOV-min are very en- 
couraging, various improvements are anticipated in the 
near future. Those include the use of a large nuclear 
database containing deformed nuclear states, advanced 
optimization protocol as in Refs. [20] HD> an d improved 
description of neutron star radius data by considering 
more realistic crust models. This work will be carried 
out under the Nuclear Low Energy Computational Ini- 
tiative (NUCLEI) [65]. 

In summary, we have successfully optimized a new 
EDF that is applicable for both finite nuclei and neutron 
stars. NS data is used to better constrain isovector inter- 
actions. As a result, the statistical extrapolation errors 
on predicted observables in neutron-rich nuclei are re- 
duced with TOV-min; hence, the functional is expected 
to yield more reliable predictions in the region of very 
neutron-rich heavy nuclei. 
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